Quantum entanglement and the self-trapping transition in polaronic systems 
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We revisit from a quantum-information perspective a classic problem of polaron theory in one 
dimension. In the context of the Holstein model we show that a simple analysis of quantum entan- 
glement between excitonic and phononic degrees of freedom allows one to effectively characterize 
both the small and large polaron regimes as well as the crossover in between. The small (large) 
polaron regime corresponds to a high (low) degree of bipartite quantum entanglement between the 
exciton and the phonon cloud that clothes the exciton. Moreover, the self-trapping transition is 
clearly displayed by a sharp drop of exciton-phonon entanglement. 



I. INTRODUCTION 

Almost three quarters of a century have passed since 
the concept of polaronic self-trapping was first conceived 
by Landau Q, Q , and yet many questions remain unan- 
swered regarding some of the most simplistic polaron 
Hamiltonians. In this paper we will approach this classic 
problem from the contemporary point of view of quan- 
tum entanglement [j| • The latter concept plays a central 
role in the burgeoning field of quantum information in 
that quantum entanglement represents the key physical 
resource at the basis of quantum information protocols 
P|. More recently, a growing amount of attention has 
been devoted to analyses of quantum entanglement in 
many-body systems undergoing quantum phase transi- 
tions y|. 

An emerging picture is that various entanglement mea- 
sures, such as the two-qubit concurrence and the entan- 
glement entropy of a subsystem, provide different, of- 
ten complementary physical insights onto the nature of 
many-body wave functions. Tools drawn from quantum 
information theory provide a deeper, previously unavail- 
able understanding of quantum correlations and their 
qualitative changes across boundaries of different phase- 
diagram regimes in these many-body systems. More- 
over, inspired by ideas from quantum-information theory, 
novel and efficient computational schemes for studying 
quantum many-body systems have been devised [||. 

The study to be carried out in this paper aims to exem- 
plify the ongoing efforts bridging the new field of quan- 
tum information theory and established areas of con- 
densed matter physics. We will show that an analysis of 
quantum entanglement between electronic (or excitonic) 
and phononic degrees of freedom in a polaronic model al- 
lows one to characterize the polaron phase diagram in a 
strikingly simple fashion. In particular, the self-trapping 
transition, i.e., a crossover between the so-called small 
and large polaron regimes, is readily captured by the be- 
havior of the linear entropy of the excitonic reduced den- 
sity matrix. We will show that the small (large) polaron 
regime corresponds to a high (low) degree of quantum 
entanglement between the exciton and its phononic envi- 
ronment. Moreover, the self-trapping transition is clearly 



displayed by a sharp drop of exciton-phonon entangle- 
ment. 

It is worthwhile to stress that, at variance with most 
of the other studies of this type, the particular form of 
quantum entanglement we are going to analyze is be- 
tween systems of distinct physical nature: the exciton, 
a finite-dimensional system, and the phonons, a bosonic 
bath. To emphasize this fact we will refer to the bipar- 
tite quantum entanglement between the exciton and its 
phonon environment in the polaron problem as hetero- 
entanglement. This situation bears resemblance to those 
in decoherence studies where a system under examina- 
tion, e.g., a qubit, is coupled with its environmental de- 
grees of freedom, which spoils the purity of the system 
state. From this point of view, the polaronic entangle- 
ment to be analyzed in this paper can be viewed as a 
measure of the decoherence of the excitonic (phononic) 
state induced by the coupling with the lattice phonons 
(excitons). 

This paper is structured as follows. In Section II, 
the Holstein model is introduced, and the correspond- 
ing phase diagram of polaronic self-trapping is presented. 
In Section III, we evaluate hetero-entanglement between 
the exciton and the phonon bath, and examine its re- 
lation to polaronic self-trapping. Discussions are given 
in Section IV, in which we show that excitonic superra- 
diance, viewed as a form of intra-exciton entanglement 
between spatially distinct excitonic modes, complement 
the exciton-phonon entanglement. 



II. SELF-TRAPPING AND PHASE DIAGRAM 
OF THE HOLSTEIN MODEL 

We first introduce a Frenkel-exciton model Hamilto- 
nian, also known as the Holstein molecular crystal model, 
which describes a lattice of two-level molecules interact- 
ing with a bath consisting of nuclear (intramolecular, in- 
termolecular, and solvent) degrees of freedom: 



H = ^fi n (q)Bt£? n 



J mn (c l )BlB n + HP h .(I) 
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Here B n (B^) are exciton annihilation (creation) oper- 
ators for the nth molecule, H ph is the bath (phonon) 
Hamiltonian, and q represents the complete set of nuclear 
coordinates. The excitonic operators satisfy the hard- 
core bosons relations (B^) 2 = B 2 = 0, [B„,B} n ] = S nm . 
It follows that for each site label n is associated to a 
two- level system i.e., a qubit. Notice that in terms of the 
Pauli matrices one has B^ = a + , and B = o~ . Despite 
this latter notation is the standard one in quantum infor- 
mation, in order to keep in line with the vast polaronic 
literature, we will stick to the B operators. Interested 
readers should not have any problem in translating the 
formulas in the Pauli matrices notation. 

Exciton-phonon interactions originate from depen- 
dence of molecular frequencies f2 n and the intermolecular 
couplings J m n on nuclear coordinates q. We adopt the 
Hamiltonian Eq. Q with the Einstein phonon Hamilto- 
nian H ph = Y^ n ^o&n&m where b\ creates a phonon of 
frequency u>q on site n, and we have one Einstein os- 
cillator per molecule. Exciton-phonon interactions enter 
through the nuclear coordinate influence on both molec- 
ular frequencies (diagonal coupling) and intermolecular 
interactions (off-diagonal coupling). Expanding Cl n (q) 
to first order in phonon coordinate q, the first term of 
Eq. Q reads "£ n Q n (q)B^B n = £ n Q„(q = 0)BlB n + 
H dza9 with the diagonal exciton-phonon coupling term 

H dia <> = ghco B\B n {bi + b n ), (2) 
n 

and g is a dimcnsionless diagonal coupling constant. Ex- 
panding J TOrl (q) to first order in phonon coordinates, 
we write the second term of Eq. as, for example, 

Em?n Jmn(ci)BlB n = E m ^„ -Wq = 0)B^B n + H od " 

with the transfer integral J m „(q = 0) = — J5 n>m ±i and 
the off-diagonal coupling term Q, H| 

H° d - = -thwoY^iBtBn+^bj+hWn+u-Snl) 
nl 

+B+B B _i(6| + bi){8 nl - 8 n -i,i)]. (3) 

The second term of © is the Hermitian conjugate of the 
first, and we have assumed nearest- neighbor coupling of 
the antisymmetric type with <f> a dimensionless param- 
eter controlling the off-diagonal coupling strength. Off- 
diagonal coupling may adopt various forms UE| other 
than the antisymmetric type 0, and can play impor- 
tant roles in electronic properties of solid. In the theory 
of high-temperature superconductivity, for example, it 
has been recently proposed , that off diagonal coupling- 
modulates the hopping integral of the Zhang- Rice singlet 
and the superexchange interaction, and is especially rel- 
evant in the low-doping regime. Eq. J5J and Eq. J2J, 
together with H ph and the zeroth-order intermolecular 
coupling term, result in the generalized Holstein Hamil- 
tonian H GH (The original Holstein Hamiltonian contains 



diagonal coupling only) [Icj : 

H GH = J2^n(ci = 0)BiB n + H dla9 

n 

+ Y. J ™^ = Q ) B ™ B n + H° d +H ph . (4) 

inn 

There are two competing energy scales in the Hol- 
stein Hamiltonian in the absence of off-diagonal exciton- 
phonon coupling, namely, the lattice relaxation energy 
g 2 uJo and the bare exciton bandwidth 4 J. Their ratio will 
be denoted the coupling strength k = g 2 u>o / 4 J which de- 
termines the size of the polaron as well as exciton-phonon 
correlations. In typical molecular crystals, g 2 < 1, in ionic 
crystals g 2 is large compared to unity, and in semicon- 
ductors g 2 is between the former two. In anthracence, for 
example, k is about 0.4, and in pyrene, about 0.8 ~ 1.6 
[Tl) . For strong exciton-phonon coupling (n 3> 1), so- 
lutions of the Holstein Hamiltonian are known as small 
polarons because the exciton-induced lattice distortion is 
confined to essentially a single exciton site 01 ■ For weak 
exciton-phonon coupling (k <C 1), the spatial extent of 
the lattice distortion is significantly increased and the re- 
sulting phonon-dressed exciton is called a large polaron. 
The crossover from a large polaron to a small polaron 
with increasing exciton-phonon coupling (essentially of- 
ten called the self-trapping transition) is rather abrupt 
for large intermolecular coupling J. In the limit of slow 
lattice motions, adiabatic polaron theories admit approx- 
imate solutions in the form of solitons. 

The one-dimensional Holstein Hamiltonian with diag- 
onal and off-dia gon al coup ling to Einstein phonons has 
been previously [lj, Ej, Ell [l6j modeled by a variational 
wave function pioneered by Toyozawa (labeled as the 
Toyozawa ansatz in Ref. 15]): 

\K)=N- 1 ^ iKn \^n)J2^-n B t\°)e ■ (5) 
n m 

Here \K) is the lowest energy polaron state with momen- 
tum K, |0) c is the exciton vacuum state, and |A^) are 
phonon wave functions centered at site n containing a co- 
herent state on each site with a displacement X^ 2 _ n : 

|A*> = CX P [- J2( X n 2 -nbl - A£*_„&n 3 )]|0) ph . (6) 

"2 

|0) p h is the phonon vacuum state, and |A^) is different 
from |Af/) only by a shift oin — n' lattice constants. The 
parameters Xf^ and ^ are obtained variationally. The 
phonon wave functions |A„ ) represent a lattice distor- 
tion forming a potential well centered at n and trapping 
the exciton with an amplitude distribution of ip, f . The 
Toyozawa Ansatz state (jHJ is not normalized: (K\K) = 
Y, nm e- iKn ^m-nSn where S* is the Debye- Waller 
factor: S% = <A£|A£_„> = exp^- 1 £ ? |Af | 2 (e^»-l)] 

and Ag , are the Fourier transform of X„ . The varia- 
tional methods are shown to be rather efficient while 
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FIG. 1: The J-g phase diagram for the Toyozawa Ansatz. A 
thin sword-shaped regime identifies with the polaronic self- 
trapping line, below which the effective mass of the polaron 
decreases drastically as g is decreased. The top (solid) line 
of the sword-shaped regime indicates the onset of bifurcation 
at K = when the polaronic structure goes through a sud- 
den change as one travels vertically downward in the phase 
diagram, and the bottom (dashed) line indicated where the 
state of the highest K (in the vicinity of K — 0) for which 
a discontinuous change in the polaronic structure is observed 
acquires two solutions as one travels vertically upward. 



remaining quantitatively accurate compared with calcu- 
lations involving far more expensive computational re- 
sources [13, IS- 

We introduce the phonon-traced exciton density ma- 



trix for the state \K): 



p« := (K\K)-'Tv ph (\K)(K\) 



(7) 



where Tr p h stands for tracing over the phonon degrees 
of freedom. To calculate we assume that exciton- 
phonon coupling leads to the formation of bands of collec- 
tive exciton-phonon states, and the many-body polaron 
wave function for the ground state band is given by the 
Toyozawa Ansatz plfllll5l |: 



(Pe)r 



,iK{n —n) I K / K* a 
l rm—nV J m'—n"- J ; 
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In this paper we confine our attention to the case of 
diagonal exciton-phonon coupling only. Localization of 
(p£~°)mm' is determined by the combined effect of elec- 
tronic confinement ip„ =0 and the overlap of the adja- 
cent nuclear wave functions S^ =0 . An efficient iterative 
procedure is employed to identify optimized Xf and i/j^ 
|14j. Computation typically starts from the J=0 axis 
where exact solutions exist. A convergent solution at one 
point in the phase diagram is used to initiate a variation 
at a neighboring point. Reversibility of the calculation 
and uniqueness of the solution are checked via altering 
initializations of the iterative procedure as the param- 
eter space is mapped. A phase diagram for the Toy- 
ozawa Ansatz spanned by J and g is presented in Figure 
1 which displays a thin sword-shaped region identified 



with the polaronic self-trapping. The top (solid) line of 
the sword-shaped regime indicates the onset of bifurca- 
tion at K — when the polaronic structure goes through 
a sudden change as one travels vertically downward in 
the phase diagram, and the bottom (dashed) line indi- 
cated where the state of the highest K (in the vicinity of 
K = 0), for which a discontinuous change in the polaron 
structure is observed, acquires two solutions as one trav- 
els vertically upward. Outside the sword-shaped area, so- 
lutions to self-consistent variational equations are unique 
and independent of how the iterative procedure is ini- 
tialized. The polaronic structure for (J, g) points above 
(below) the sword-shaped area is traditionally identified 
as small (large) polarons. 



III. HETERO-ENTANGLEMENT AND 
POLARONIC SELF-TRAPPING 

For a finite bare exciton band width J > 0, the lowest 
energy state for diagonal coupling only has zero crystal 
momentum K = 0. In this lowest-energy K = state, 
hetero-entanglement between the two species in the Hol- 
stcin Hamiltonian, the exciton and the phonons, as mea- 
sured by the linear entropy, has the form 



B:=l-Tr e [(pf=Y] 



(8) 



The function JSJ is a linearized version of the von Neu- 
mann entropy S(p) = — Trpln p; it shares with this latter 
quantity the properties: (i) E = p = p , i.e., p is 
a pure state; (ii) -E max = 1 — l/D = E(t/D), i.e., the 
linear entropy is at its maximum for the totally mixed 
state l/D (D = is the dimension of the space). More- 
over, for qubits, the linear and von Neumann entropy are 
monotonic functions of each other. The linear entropy 
JSJ) has a close relation with the so-called 2-Renyi en- 
tropy 0| . These considerations show that (JBJ) represent 
a good choice as a simple entanglement measure. Cal- 
culated hetero-entanglement is shown in Figure [21 for the 
entire J vs g phase diagram. It is important to note that 
the qualitative behavior of the ground state entanglement 
as a function of control parameters discussed above does 
not depend on the specific choice of the linear entropy as 
the entanglement measure; analogous calculations per- 
formed with the von Neumann entropy give rise to the 
very same qualitative picture [2(j 

On the J = vertical axis, the exact wave 
function of the Holstein polaron can be written as 
E„S, t jO) c |AJ=°),where (A^A^ ) = 5 mn . Therefore, 
the hetero-entanglement between the two species reaches 



its maximum on the J = vertical axis: E n 



1-N- 







On the horizontal axis g — 0, on the contrary, the K 
state is separable with respect to the two species. 

As shown in the upper panel of Figure [21 exciton- 
phonon hetero-entanglement forms a cliff. The edge of 
the cliff is found overlapping with the thin sword-shaped 
area in Figure 1 (i.e., the self-trapping line), and can 
be fitted empirically by J c — (g c — l) 2 hu>o: on the cliff 
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linear entropy of the excitonic (or phononic) reduced den- 
sity matrix. The latter has been explicitly computed by 
resorting to utilization of the Toyozawa Ansatz. Diffi- 
culties associated to the infinite dimensionality of the 
phonon degrees of freedom can be circumvented by us- 
ing the Toyozawa variational Ansatz to approximate the 
ground-state polaronic wave function, for which an exact 
form remains elusive. This choice of Ansatz states also 
allows us to carry out a large portion of the calculations 
in an analytical, conceptually transparent fashion. The 
study of the entanglement behavior, as a function of the 
controlling parameters J and g (exciton hopping ampli- 
tude and exciton-phonon coupling respectively), allows 
a very simple characterization of the zero-temperature 
phase diagram of the model. The self-trapping transi- 
tion from large to small polaron regimes can be under- 
stood in terms of a sharp increase of the exciton-phonon 
entanglement. 

We also note that this form of hetero-entanglement 
is complementary to the superradiant behavior of the 
excitonic system. Superradiance (coherent spontaneous 
emission) is the enhanced radiative decay compared to 
that of a monomer as a result of the coherent nature 
of the electronic excited states . If all transition 

dipoles of the A monomers are parallel, superradiance of 
the K = state can be calculated from 
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FIG. 2: Hetero-entanglement between the exciton and the 
phonons as measured by the linear entropy 1 — Tr o [(p^ =0 ) 2 ] is 
displayed in the upper panel for the entire J-g phase diagram. 
The entanglement is calculated for the lowest energy state 
with zero crystal momentum K=0. The solid line in the lower 
panel is the edge of the cliff which separates the large and 
small entanglement regions. 



plateau, E approaches 1 — A -1 ; on the other side of the 
cliff edge, E decreases rapidly to zero with increasing J 
or decreasing g. In the lower panel, the solid line is the 
projection of the cliff edge onto the J-g plane which sep- 
arates the large and small entanglement phases. The 
usual distinction between the small and large polarons 
can thus be rephrased as follows: the small polaron is 
a maximally entangled exciton-phonon entity, while the 
large polaron has much-reduced exciton-phonon hetero- 
entanglement. The hetero-entanglement is therefore a 
good measure of large and small polarons and the tran- 
sition in between. 



IV. DISCUSSIONS 

In this paper we have revisited the classic polaron 
problem utilizing the concept of quantum entanglement. 
The exciton-phonon coupling, described by the Holstein 
Hamiltonian, induces strong quantum coherences be- 
tween the two heterogeneous physical degrees of freedom. 
These coherences can be conveniently quantified by the 
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(9) 



It follows that excitonic superradiance and the exciton- 
phonon hetero-entanglement are complementary. For ex- 
ample, on the vertical J=0 axis in Figure 2, the linear 
hetero-entanglement JSJ is at its maximum, 1 — A -1 , 
while the corresponding K = superradiance reaches its 
minimum value 1; on the horizontal g — axis, the linear 
hetero-entanglement vanishes while corresponding K = 
superradiance reaches its maximum value A. For a given 
transfer integral J, as the exciton-phonon coupling g 
is reduced, the superradiance gains while the exciton- 
phonon hetero-entanglement decreases. Superradiance 
can be regarded as a form of the quantum coherence 
between different components of the exciton wave func- 
tion localized on different sites. This coherence, in turn, 
can be viewed as a form of entanglement between spa- 
tially distinct excitonic modes (for a definition of mode 
entanglement, see, e.g., Hj|). Thus, the aforementioned 
complementarity phenomenon can be viewed as a form 
of entanglement transfer from the inter-exciton-phonon 
type to the intra-exciton mode entanglement, and vice 
versa. 

The role of off-diagonal exciton-phonon coupling as 
well as analyses of entanglement of polaronic states with 
nonzero crystal momenta are the subject of ongoing in- 
vestigations. The Toyozawa Ansatz is again well-suited 
to capture the underlying physics in the presence of off- 
diagonal exciton-phonon coupling. We believe that the 
results presented in this paper exemplify the fact that 
implementations of notions and apparatuses drawn from 
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the emerging field of quantum information science are in- 
valuable for gaining fresh insights for classic problems in 
condensed matter physics. 
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